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ABSTRACT 

Context. Accretion and ejection are tightly connected and represent the fundamental mechanisms regulating star formation. However, 
the exact physical processes involved are not yet fully understood. 

Aims. We present high angular and spectral resolution observations of the Bry emitting region in the Herbig Ae star HD 163296 
(MWC 275) in order to probe the origin of this line and constrain the physical processes taking place at sub-AU scales in the circum- 
stellar region. 

Methods. By means of VLTI-AMBER observations at high spectral resolution (R~12000), we studied interferometric visibilities, 
wavelength-differential phases, and closure phases across the Bry line of HD 163296. To constrain the physical origin of the Bry line 
in Herbig Ae stars, all the interferometric observables were compared with the predictions of a line radiative transfer disc wind model. 
Results. The measured visibilities clearly increase within the Bry line, indicating that the Bry emitting region is more compact than 
the continuum. By fitting a geometric Gaussian model to the continuum-corrected Bry visibilities, we derived a compact radius of the 
Bry emitting region of ~0.07±0.02 AU (Gaussian half width at half maximum; or a ring-fit radius of ~0.08±0.02 AU). To interpret the 
observations, we developed a magneto-centrifugally driven disc wind model. Our best disc wind model is able to reproduce, within the 
errors, all the interferometric observables and it predicts a launching region with an outer radius of ~0.04 AU. However, the intensity 
distribution of the entire disc wind emitting region extends up to —0.16 AU. 

Conclusions. Our observations, along with a detailed modelling of the Bry emitting region, suggest that most of the Bry emission in 
HD 163296 originates from a disc wind with a launching region that is over five times more compact than previous estimates of the 
continuum dust rim radius. 

Key words, stars: formation - stars: circumstellar matter - ISM: jets and outflows - ISM: individual objects: MWC275, HD163296 
- Infrared: ISM - techniques: interferometric 


1. Introduction 


As a general view, the accretion and ejection processes in young 
stellar objects (YSOs) proceed as part of the matter in the disc 
is lifted up along the magnetospheric accretion columns and 
accretes onto the stellar surface, whereas a small amount of 
the accreting matter is instead accelerated and collimated out¬ 
wards in the form of protostellar winds or jets. Accretion and 
ejection are tightly connected and represent the fundamental 
mechanisms regulating the formation of a star. The first indi¬ 
cations of such a connection were found in Classical T Tauri 
stars (CTTSs) and Herbig AeBe stars through the discovery of 
a correlation between the luminosity of jet line tracers (such 


as the fOil 6300 A l ine) and the infrared excess luminosity 
(lHartigan et alJfl995L ICorcoran & Ravlfl998l) . The exact physi¬ 
cal processes driving and connecting the accretion and ejection 
activity are, however, not fully understood. One of the main rea¬ 
sons is the complex structure of the accretion-ejection region 
and the small angular scales involved: within less than 1 AU 


Send offprint requests to: rgarcia@cp.dias.ie 
* Based on observations collected at the European Southern Obser¬ 
vatory Paranal, Chile (ESO programme 089.C-0443(A)). 


from the central source (i.e. <10mas at a distance of lOOpc) 
emission from the hot gaseous disc, outflowing material, and the 
accretion columns is expected. In this context, near-infrared in¬ 
terferometry and, in particular, HI Bry spectro-interferometric 
observations, have become an essential tool with which to in¬ 
vestigate the hot and dense gas at sub-AU scales. The first 
Bry spectro-interferometric observations of Herbig AeBe stars 
showed that in some sources the bulk of the Bry emission arises 
from an exte nded component (probably tracing a disc or stel- 
lar wi nd; e.g. Malbet et td1l2007l: lEisner et aPl2007l: iKraus et aD 
[200& 1 Weigelt et al.120 111) . in contrast with the scenario in which 
the Bry line would be related with magnetospheric accretion, 
and so unresolved by current IR interferometers. Despite its 
unclear origin, the Bry line luminosity is directly correlated 
with the accretion luminosity through e mpirical relationships 
dMuzerolle et alJfl998t ICalvet et al.ll2004l) and extensively used 
to estimate the mass accretion rates of YSOs with different 
masses (0.1 -10Mo) and at different evolutionary stages (10 5 - 
10 7 yr; e.g. Garcia Lopez et al.l 120061: lAntoniucci et alJ 120111: 
ICaratti o Garatti et al.ll2012). 

In order to further constrain the origin of the Bry line in 
Herbig AeBe stars and probe the accretion-ejection connec- 
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tion, we have started a detailed high spectral resolution (HR) 
VLTI/AMBER study of Herbig AeBe stars. We present our re¬ 
sults on the Herbig Ae star HD 163296 (MWC275). This star is 
one of the best candidates with which to perform a high angular 
and spectral study of the accretion-ejection region: HD163296 
is located at just ~ 119 pc (Ivan Leeuwenl 120071) . it is among 
the few Herbig Ae stars driving_a well collimated bipolar jet 
(HH409; see e.g. WasseH_gtaLll20Q6 ), and it is very likely a sin¬ 
gle star (ISwartz et al. 2005 : Montesinos et alJ[2009l) . The main 
HD 163296 stellar parameters are listed in TableQ] 

The protoplanetary disc around HD 1 63296 has been ex¬ 
tensively stu died (e.g. Ilsella et al.l 20071 : iBenistv et al.l 120101: 
iTilling et alJl2012t iRosenfeld et al. 1 201 3l) . Recently, ALMA ob¬ 
servations have revealed a vertical temperature gradient across 
the disc, suggesting dust settling and/or migration (se e e.g. 
Ide Gregorio-Monsalvo et all 1201 3c iMathews et al.l l2013h . De¬ 
spite the evidence of dust growth, the disc still has enough gas 
content to drive a CO disc w i nd and a large-scale collimated 
bipolar jet (IWassell et alJl2006LlKlaassen et al.ll2013l) . 

Previous, spectro-interferometric studies of HD 163296 suc¬ 
cessfully resolved the inner edge of the circumstellar disc, 
locating it at distance valu es from the star ranging from 
~0.19AU to ~0.45AU ( e.g. RenardetaL 20101; Jenistv et al. 
l2010t lEisner et~akl l2010l 1201 4t iTannirkulam et al. 2008allbl) .~ft 
is withi n this region that the extended Bry emission is being 
emitted dKraus et alJl2008t lEisner et al .1120141) . although, the ex¬ 
act origin of the line emission, and the nature of the innermost 
gaseous disc are still unclear. 

In the following, we will present our VLTI/AMBER HR ob¬ 
servations (R~ 12 000) across the Bry line of HD 163296. The 
observations, data reduction, and calibration, as well as the in¬ 
terferometric observables derived from our observations are pre¬ 
sented in Sect. [2] and Sect. [3] In Sects. [4] [5] and [6] a description 
of the models used to interpret our results and their comparison 
with the observations are shown. Finally, in Sect. [8] a summary 
of our results and the main conclusions are presented. 


2. Interferometric VLTI/AMBER/FINITO observations 
with spectral resolution of 12000 

We observed HD 163296 during three nights in May and June 
2012 with the ESO Very Large Telescope Interferometer (VLTI) 
and its AMBER beam combiner instrument dPetrov et alil2007l) . 
For these observations, we used the three 8 m Unit Telescopes 
UT2, UT3, and UT4 with AMBER in the high spectral resolution 
mode (R=12 000). The observation parameters are described in 
Table[2] To obtain interferograms with a high signal-to-noise ra¬ 
tio (S/N), we used the fringe tracker FINITO for cophasing and 
a detector integration time of 6.0 s per interferogram. 

For data reduction, we used our own data reduction soft¬ 
ware based on the P2VM algorithm dTatulli et al . 12007 1 to de¬ 
rive wavelength-dependent visibilities, wavelength-differential 
phases, and closure phases. Along with the science observations 
of HD 163296, we observed the interferometric calibrator stars 
HD 163955, HD 162255, and HD 160915, which were used to 
calibrate the transfer function. The fringe tracking performance 
of FINITO is usually different during the target and calibra¬ 
tor observations. Therefore, we used all published archival low- 
resolution AMBER observations of HD 163296 to calibrate the 
continuum visibilities. The errors of the K-band squared visibi¬ 
lities are approximately +5% or slightly better during the best 
nights. Each of our three recorded datasets consists of 75 tar¬ 
get and calibrator interferograms. Because all three observations 


velocity (LSR) [km/s] 



2.160 2.162 2.164 2.166 2.168 2.170 2.172 

X [pm] 


Fig. 1. AMBER observation of HD 163296 with three different base¬ 
lines (45.0 m: red upper line, 57.5 m: green middle line, and 81.6 m: 
blue lower line) and spectral resolution R=12 000. From top to bottom: 
wavelength dependence of flux, visibilities, wavelength-differential 
phases, and closure phases. For clarity, the differential phases of the 
first and last baselines are shifted by +30 and -30°, respectively. The 
results shown correspond to the averaged value over three different ob¬ 
servations (see Sect.[2]and Appendix for details). The wavelength scale 
is with respect to the local standard of rest (LSR) and corrected by a 
cloud velocity of ~8 km s’ 1 dlsella et aDl2007h . 


were performed with very similar projected baseline lengths and 
position angles (PAs), we averaged the results of all three obser¬ 
vations in order to get a better S/N per spectral channel (Fig.Q}. 
The smaller errors of the averages of the three datasets were ob¬ 
tained from the errors of the individual results in the standard 
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Table 1. HD 163296 stellar parameters 


SpT (1) 

r r , (2) 

eff 

d< 3 > 

r1 2) 

M?> 

m <4) 

LVX acc 

M (b) 

iV1 out 


(K) 

(pc) 

(Ro) 

(Mo) 

(lO- 7 M 0 /t/r) 

(M Jyr) 

Al V 

9250 

119+11 

2.3 

2.2 

0.8-4.5* 

5xl0~ 10 -2xl0~ 7 ** 


References. (11 Ivan den A ncker et all dl998l); (2 1 iMendigut ia et all 1 20131): (31 Ivan Leeuwenl (120071) : (4) iGarcia Lopez et alj ([2006); 
iDonehew & Brittainl d201 lh i lMendigutfa et al.l 120 1 3l) : Gt lEllerhroek et ail 2014f) : lKlaassen et alJ(l2013h 

* Range of M acc values reported in the literature (see references in table) derived from the luminosity of the Bry line. 

** The first and second values correspond to the atomic ([S 11 ], [0 1 ]) and molecular (CO) jet/outflow components. 


way. Because the absolute continuum visibilities were calibrated 
using the published archive low resolution (LR) observations, 
in the next step, we took the errors of these LR archival data 
into account and we calculated the total error in the standard 
way. The three individual results are shown in the Appendix 
(Fig. [AT]). The wavelength calibration was accomplished using 
the numerousjelluric lines present in the region 2.15-2.19pm 
(see lWeigelt et al.ll201 lL for more details on the wavelength ca¬ 
libration method). We estimate an uncertainty in the wavelength 
calibration of ~ 3 kms -1 . 


3. Interferometric observables 

Our VLTI/AMBER observations provide four direct observ¬ 
ables: Bry line profile, visibilities, differential phases, and clo¬ 
sure phases. These observables allow us to retrieve information 
about the size and kinematics of the Bry emitting region. 

FigureQ] shows (from top to bottom) the Bry line profile, 
wavelength visibilities, differential phases, and closure phase 
of our Bry spectro-interferometric observations of HD 163296. 
The wavelength-dependent visibilities (second panel from top) 
clearly increase within the Bry line in all our baselines. This in¬ 
dicates that the Bry emitting region is more compact than the 
continuum emitting region. 

Previous spectro-interferometric studies of HD 163296 at 
medium resolution had a lso detected an increase of the visibi - 
lity within the Bry line dKraus et al.l 12008b lEisner et aDl2014l) . 
Our HR mode AMBER observations have, however, allowed 
us to measure the visibilities and phases in ~30 different spec¬ 
tral channels across the Doppler-broadened Bry line, helping us 
to better distinguish the contribution of the line to that of the 
continuum. The measured differential and closure phases do not 
clearly show any significant deviation from zero within the error 
bars of ±5°, and ±15°, respectively (Fig.Q] bottom panel). 


4. Geometric modelling: The characteristic size of 
the Bry line-emitting region 

The high spectral resolution of the AMBER data allows us to 
measure the visibilities of the pure Bry line-emitting region 
for several spectral channels across the Bry emission line (see 
Fig.ED. These continuum-compensated line visibilities (re¬ 
quired for the size determination of the Bry line-emitting re¬ 
gion) were calculated in the following way. Within the wave¬ 
length region of the Bry line emission, the measured visibi¬ 
lity has two constituents: the pure line-emitting component and 
the continuum-emitting component, which includes continuum 
emission from both the circumstellar environment and the un¬ 
resolved central star. The emission line visibility Vury in each 


spectral channel can be written as 


F BrykBry - 

- Vi^ 


'totVtotP + IFeVcl 

= |F tot V tot - F C V C |, 


■ 2 FtotVtot FrVr • CO id) 


( 1 ) 
( 2 ) 

if the differential phase <p is zero (IWeigelt et al.ll2007l) . The value 
F Bry is the wavelength-dependent line flux, V t ot (F tot ) denotes the 
measured total visibility (flux) in the Bry line, V c (F c ) is the visi¬ 
bility (flux) in the continuum, and <!> is the measured wavelength- 
differential phase within the Bry line. The intrinsic Bry photo- 
spheric absorption feature was taken into account in this ana¬ 
lysis by considering a synthetic spectrum with the same spec¬ 
tral type and surface g ravity of HD 163296 (lEisner et alJl2010t 
iMendigutfa et al.ll2013l) . 

We have derived the pure line visibility only in spectral re¬ 
gions where the line flux is higher than ~10% of the continuum 
flux. The results are shown in Fig. iBTl The average line visibi¬ 
lity is ~ 1 at the shortest baselines (45 m), ~0.9 at medium base¬ 
lines (57.5 m), and ~0.8 at the longest baselines (81.6 m). The 
approximate size of the line-emitting region was obtained by fit¬ 
ting a circular symmetric Gaussian model to the line visibilities 
presented in Fig. IBTl We obtained a Gaussian half width at half 
maximum (HWHM) radius of 0.6+0.2 mas or ~0.07+0.02 AU. 
Similarly, by fitting a ring model with a ring width of 20% 
of the inner ring radius, an inner radius of 0.7+0.2 mas or 
~0.08+0.02 AU, is found. This radius is smaller than the in¬ 
ner continuum dust rim radius of R„- m ~0 .1 9-0.45 AU reported 
in the literature (Tannirkulam et al. 2008a b: lBenistv et al ,ll20 1 (j 


lEisner et Id1l20l4 


5. Disc and wind model 

Following our previous work on the Herbig Be star MWC297, 
we developed a disc wind model to explain the observed 
Bry line pr ofile, line visibilities , differential phases, and clo¬ 
sure phase dWehge lLetal] 1201 it iGrinin & Tambovtseval 1201 lfc 
Tambovts eva et al.l2014l) . A detailed description of the disc wind 
model, the complete algorithm for the model computation, and 
a detailed description on how the disc wind model solution de¬ 
pends on the kinematic and physical parameters can be found 
in these papers. In the following, we will briefly summarise the 
main characteristics of our disc wind model. 

The disc wind model employed here is based on the 
magne to-centrifugal disc wind model of iBlandford & Pavnel 
( 1982 ;). adapted for the particular case of Herbig AeBe stars. For 
simplicity, the disc wind consists only of hydrogen atoms with 
constant temperature (~10000K) along the wind streamlines. 
This approximatio n is in agreement with the so -called warm 
disc wind models (ISafierll 19931: IGarcia et al.ll2001l) . in which the 
wind is rapidly heated by ambipolar diffusion to a temperature 
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Table 2. Log of the VLTI/AMBER/FINITO observations of HD 163296 and its calibrators. 


HD 163296 
Observation 
Date 

Time [UT] 
Start End 

Unit Telescope 
array 

Spectral 

mode__ 

Wavelength 

range 

N 

DIT" 

[s] 

N_ 

Seeing 

[arcsec] 

Calibrator 

index 

(see below) 

2012 May 11 

08:41 

08:49 

UT2-UT3-UT4 

HR-K-F 

2.147-2.194 

6 

75 

0.6-0.8 

(a)+(b) 

2012 May 12 

09:20 

09:28 

UT2-UT3-UT4 

HR-K-F 

2.147-2.194 

6 

75 

0.9-1.5 

(c) 

2012 fun. 04 

07:23 

07:30 

UT2-UT3-UT4 

HR-K-F 

2.147-2.194 

6 

75 

0.7-0.9 

(d)+(e) 


Calibrator 

Name 

(index) 

Date 

Time [UT] 
Start End 

UT array 

N“ 

Seeing 

[arcsec] 

Uniform-disc 

diameter^ 

[mas] 

HD 162255 (a) 

2012 May 11 

08:09 

08:17 

UT2-UT3-UT4 

75 

0.6-0.7 

0.34+0.02 

HD 160915 (b) 

2012 May 11 

09:05 

09:13 

UT2-UT3-UT4 

75 

0.6-0.7 

0.68+0.05 

HD 160915 (c) 

2012 May 12 

08:56 

09:04 

UT2-UT3-UT4 

75 

1.0-1.5 

0.68+0.05 

HD 162255 (d) 

2012 Jun. 04 

06:35 

06:42 

UT2-UT3-UT4 

75 

0.9-1.2 

0.34+0.02 

HD 163955 (e) 

2012 Jun. 04 

07:45 

07:53 

UT2-UT3-UT4 

75 

0.6-0.8 

0.40+0.03 


Note: Object and calibrators were observed with the same spectral mode, DIT, and wavelength range. 
a Detector integration time per interferogram. 
b Number of interferograms. 

c High spectral resolution mode in the K band using the fringe tracker FINITO. 
d UD diameter taken from JMMC Stellar Diameters Catalogue - JSDC dLafrasse et al.l20loh . 


of ~10000K. In these models, the wind electron temperature in 
the acceleration zone near the disc surface is not high enough to 
excite the Bry line emission. Therefore, in our model, the low- 
temperature region below a certain height value does not con¬ 
tribute to the Bry disc wind emission. 

For the calculations of the model images of the emitting re¬ 
gion in the line frequencies and their corresponding interfero¬ 
metric quantities (i.e. visibilities and phases), we used the same 
approach as in IWeigelt et al.l (1201 ll) . For the calculations of the 
ionization state and the number densities of the atomic levels, we 
adopte d th e numerical codes d e velope d bv lGrinin & Mitskevichl 
(Il990h and iTambovtseva et all (1200 ll) for movin g media. These 
codes are based on the Sobolev approximation (tSobolevIl 19601) 
in combination with the e xact integration of the line intensities 
(see Appendix A in lWeigelt et akKoTTl) . This method takes into 
account the radiative coupling in the local environment of each 
point caused by multiple scattering. The calculations were made 
for the 15-level model of the hydrogen atom plus continuum, 
taking into account both colli sion and radiative proces ses of ex¬ 
citation and ionization (see lGrinin & MitskevichlfT990l for more 
details). 

Very briefly, our disc wind model considers a wind launched 
from an inner radius u>\ to an outer radius w,y (wind foot- 
points). The disc wind half opening angle (6) is defined as the 
angle between the innermost wind streamline and the system 
axi s. The angle range ( 30°-45 °) in Tabled is consistent with 
the iBlandford & Payn el d 1982b disc wind solution. The local 
mass-loss rate per unit area on the disc surface is defined as 
m(to) ~ u>~ J , where y is the so-called mass-loading parameter 
that controls the ejection efficiency. Thus, the total mass-loss rate 
is given by 

(i)N 

M w — 2 J' m w (u>) 2 7r u> du>. (3) 

In all our disc wind models, the accretion disc is assumed 
to be transparent to radiation at radii up to the dust sublimation 


radius. At larger radii, the disc is an opaque screen that shields a 
large fraction of the wind emission because the observer cannot 
see the disc wind emission on the back side of the disc. It should 
be noted that this assumption (i.e. opaque disc with an inner gap) 
is only valid if the mass accretion rate ( M nrr ) does not exceed a 
value of ~1(T 6 M Q yr~' dMuzerolle et al.ll2004l) . 

All free model parameters are listed in Tabled Some model 
parameters such as M w , the stellar parameters and the disc in¬ 
clination angle were taken from the literature (see Tables[3] and 
Q}. The modelled total mass-loss rate (M„F1) was allowed to vary 
within ~0.1-1.0 times the M acc observed average value (see Ta¬ 
ble!]]). In addition, the estimated value of the corotation radius is 
set as a lower limit to the inner disc wind launching radius value 
u)\ (see Sect. 01. These parameters, along with our interferome¬ 
tric observations further constrain our disc wind model solution. 
For instance, the measured visibilities constrain the size of the 
disc wind launching region, and the adopted M w and measured 
line p rofile constrain the wind dynamics (see ITambovtseva et al.l 
12014 for a complete discussion). 

Finally, in order to compute the model differential phases, 
a continuum model consisting of tw o components, an inner 
and an outer disc, was adopted (e.g. iTannirkulam et alJ l2008a: 
iBenistv et al.l20 10). This model is able to reproduce both the ob¬ 
served continuum visibilities and continuum flux of HD 163296. 
The dominant outer disc is a temperature-gradient disc with a 
power law index of -0.50, inner temperature of 1500 K, and a 
dust rim radius of 20 R* (i.e. ~0.21 AU). The stel lar flux is 0.14 
of the observed flux in the K-band dBenistv et al.ll2010h . The ad¬ 
ditional optically thin inner disc has an inner radius of 5 R*, an 
outer radius equal to the inner radius of the outer disc (i.e. 20 R»), 
and a constant intensity distribution flux of 20% of the observed 
continuum flux. 


1 From now on, we use M w for the modelled total mass-loss rate and 
M 0 u t for the observed value from the literature. 
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Table 3. Disc wind model parameters 


Parameters^ 

RangU 

MW6* 

MW26 

MW47 

Temperature (K) 

8000-10000 

10000 

10000 

10000 

Half opening angle (6) 

30° - 45° 

45° 

30° 

20° 

Inner radius (cuBR*)) 

2-3 

2.0 (0.02 AU) 

2.0 (0.02 AU) 

5.0 (0.05 AU) 

Outer radius (a>jv(R*)) 

4-30 

4.0 (0.04 AU) 

15.0(0.16 AU) 

15.0 (0.16 AU) 

Acceleration parameter (J3 ) 

1-7 

5 

3 

4 

Mass load parameter (y) 

1-5 

3 

3 

3 

Mass loss rate (M„,(M 0 /yr)) 

1 

O 

1 

oo 

1 

O 

5xl0~ 8 

5xl0~ 8 

5xl0^ 8 


a The definition of the listed disc wind model parameters is described in detail in lWeigelt et alj d2011h . 

b Disc_inclination and position angle (major axis) of ~38° and ~138°, respectively, were assumed (gaseous ALMA disc; 

Ide Gregorio-Monsalvo et alj|2013i) . as well as the stellar parameters listed in TableQ] 


6. Results 

In order to find a disc wind model that is able to approximately 
reproduce all interferometric observables, we calculated more 
than 100 models within the parameter ranges listed in Tabled 
All our models were computed assuming a disc inclination (an¬ 
gle between the line of sight and the system axis of the gaseous 
ALMA gas disc) of ~38° and a posi tion angle of the major 
axis of the disc of ~138° (Ide Gregorio-Monsalvo et ahll2013l) . 
We computed the intensity distributions of the models for all the 
wavelength channels across the Bry line. 

Because our observations do not provide us with wavelength- 
dependent images (as the modelling does), but only with the Bry 
line profile and information on visibilities and phases, we first 
had to derive the following interferometric quantities from the 
wavelength-dependent model intensity distributions: model line 
profiles, model visibilities, and differential and closure phases of 
the model. All the interferometric model quantities were com¬ 
puted for the same baseline lengths, baseline PAs, and spectral 
resolution as the observations. This method allowed us to com¬ 
pare the observed interferometric quantities (i.e. line profile, vi¬ 
sibilities, differential phases, and closure phases) with the corres¬ 
ponding interferometric model quantities and to find out which 
of our models can best reproduce the observations. 

Figures[2] and IB. II show a comparison of the observations 
with the interferometric quantities derived from our best disc 
wind model MW6 (see Tabled for a list of the model parame¬ 
ters). The intensity distribution maps of this model at several 
radial velocities are presented in Fig. [3] All model parameters 
are listed in Tabled As shown in the figures, model MW6 ap¬ 
proximately agrees with the observed line profile (Fig. El up¬ 
per panel), visibilities (Fig. [2] middle panel), and small differ¬ 
ential phases and closure phases (Fig. El the measured differen¬ 
tial and closure phases are approximately zero within our errors; 
the model closure phase is ~1°). Moreover, this model is also 
able to approximately reproduce the visibilities and small dif¬ 
ferential ph ase of previous Keck inte rferometric observations of 
this source (lEisner et al.ll201C)L 120141) . However, it slightly over¬ 
predicts the pure Bry line visibilities for the longest baseline 
(see Fig. IB.11 1. This can also be seen in Fig. El where the model 
continuum visibility is slightly below the observed value for the 
longest baseline. 

In order to obtain a lower value of the modelled pure line 
visibility at the longest baseline, it would be necessary to in¬ 
crease the size of the disc wind region (e.g. lon). Two examples 
of such models (i.e. with ojn >4R„) can be found in Fig. PI 
(model MW26) and Fig. O (model MW47). Model MW26 has 
the same o>\ value as our best model MW6, but a u>n value 
almost 4 times larger (0.16 AU). Model MW47 has larger val¬ 


ues of a>i (0.05AU) and u>^ (0.16AU). The modelled interfe¬ 
rometric results presented in Figs. PI and Figs. 1C.21 show that 
these larger disc wind launching radii have little influence on the 
line visibility. Both new models have, however, large differential 
phases (>34°) that disagree with the observed differential phases 
(approximately zero within the average error of ~ +5°). These 
large differential phases are caused by the larger radii, that pro¬ 
duce larger photo centre shifts with respect to the continuum. 
It should be noted that the acceleration parameters of models 
MW6, MW27, and MW47, as well as the half opening angles, 
differ slightly from model to model (/?=5, 3, and 4, and 0=45°, 
30°, and 20°, respectively). However, these model p aram eters 
mainl y modify the Bry line profile (see e.g. Tambovtsev a et al.l 
120141) . and have little effect on the line visibilities and differ¬ 
ential phases which mostly depend on the a>\ and w,y parame¬ 
ters. Therefore, a good match between our interferometric ob¬ 
servables (line profile, visibilities, and differential and closure 
phases) and the modelled results can only be achieved when 
combining lower values of the inner and outer disc wind foot- 
points, that is, when a compact disc wind is considered. 


7. Discussion 

The disc wind launching region of our best model MW6 is quite 
compact, extending from an inner radius of ~0.02 AU to an outer 
radius of ~0.04 AU. This latter value is not far from the lowes_t 
outer launching radius of 0.03 AU derived from iFerreira et al.l 
(120061) from a comparison of optical jet observations of TTauri 
stars (typical R co ,. ~0.1 AU) and the predictions from extended 
disc wind models. Even if a direct comparison between TTauri 
and Herbig AeBe stars is difficult, this result suggests that disc 
wind models with low values can explain some observations 
of extended emission in YSOs. However, Fig. [3] shows that the 
entire disc wind emitting region extends to an outer radius of 
~15R« or ~0.16AU. As discussed in the previous section, this 
allows us to obtain approximately the required low visibilities at 
the longest baseline, as well as small differential phases, and to 
fit the Bry line profile. 

It should be noted that our modelling does not exclude 
that a fraction of the Bry flux might be emitted from a com- 
pact magnetosphere, a st ellar wind, and/or an X-wind (see e.g. 
Tambovtseva et al. 2014, for examples of hybrid models). Be¬ 
cause of the small magnetic fields measured in Herbig AeBe 
stars, the distance at which the circumstellar disc of Herbig AeBe 
stars is truncated by the stellar magnetic field (the truncation 
radius) is expected to be much smaller than in CTTSs and it 
is located within the corotation radius (R cor ) (IShu et al.l Il994t 
iMuzerolle et al.ll2004l) . For the case of HD 163296, the corota- 
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Fig. 2. Comparison of our interferometric observations with the interfe¬ 
rometric quantities derived from our best disc wind model MW6. From 
top to bottom, observed Bry line profile (grey) and model line profile 
(pink), observed visibilities (grey dots with error bars) and model vi¬ 
sibilities (coloured lines), observed and model wavelength-differential 
phases, and, observed and model closured phases. 

tion radius, R for , is ~1.6R* or -0.02 AL0. Therefore, the magne¬ 
tosphere outer radius (<0.02 AU) is much smaller than the outer 
radius of the whole disc wind emitting region (—0.16 AU; see 
Fig. [3]), and comparable to the inner launching radius parame- 


2 The corotation value of 1.6 R, was derived from the stellar pa- 
ramet ers pre s ented in Tablet]] assuming an inclination angle of 38° 
dde Gregorio-Monsalvo et aT1l2013l . see TableO . and a projected rota¬ 
tional velocity of v sin i = 133 km s _1 iMontesinos et al.ll2009h . 


ter u >i in our disc wind model MW6. As a consequence, any 
additional emission in our model located within the corotation 
radius (-0.02 AU) would increase the modelled Bry visibilities 
of the disc wind intensity distribution. To illustrate this effect, 
we employed a hybrid model (disc wind plus magnetosphere) 
as described in iTambovtseva et all (120141) . and add the contri¬ 
bution of a compact magnetosphere accounting for 40% of the 
total Bry flux to our best model MW6 (a brief description of 
the magnestosphere parameters and their values can be found in 
AnnendixlDl). To account for the extra emission from the magne¬ 
tosphere and obtain a good fit to the line profile, the Bry emitting 
region in model MW6 was reduced by increasing the size of the 
low-temperature region at the wind base where the gas does not 
produce Bry emission (see Sect. [5}. In this way, the disc wind 
Bry emission was reduced to a 60% of the total Bry flux. Fig- 
ure lD. II shows the result of this model overplotted on our obser¬ 
vations: the pure Bry line visibilities are now ~1 at all baselines. 
To conclude, our hybrid model (i.e. disc wind plus magnetos¬ 
phere) suggests that the observed low visibility values can only 
be achieved with larger disc wind footpoint radii than in model 
MW6, but this would in turn produce excessive values of the 
differential phase as well, contradicting our observations. 

Nevertheless, we would like to point out that the aim of this 
paper is to investigate which mechanism is responsible for most 
of the Bry emission in HD 163296. Therefore, our observations 
do not rule out that some of the Bry emission is produced in an 
X-wind and/or a magnetosphere in HD 163296, but they would 
probably have to act in combination with a more extended emis¬ 
sion like the disc wind reported here in order to reproduce the 
observed Bry line visibilities and the line profile. 

The presence of a disc wind in HD 163296 is also supported 
by ALMA observations which revealed an extended CO rota¬ 
ting disc wind emerging from this source (Klaas sen et alJl20l3h . 
Our observations also indicate a compact disc wind component 
traced by the Bry line, and launched well within the dust sub¬ 
limation radius. This, in turn, would explain why no dust has 
been detected along the long-scale optical/infrared bipolar jet 
(lEllerbroek et al.ll2014l) . that is, the jet is launched from an in¬ 
ner disc region where none or only little dust is present. 

8. Summary and conclusions 

In this paper, we present AMBER high spectral resolu¬ 
tion (R=12 000) observations of the young Herbig Ae star 
HD 163296. The high spectral resolution has allowed us to sam¬ 
ple the interferometric observables (i.e. line profile, visibilities, 
differential phases, and closure phases) in many spectral chan¬ 
nels across the Bry line. Our observations show that the vi¬ 
sibility within the Bry line is higher than in the continuum 
(FigJD. Therefore, the Bry emitting region is less extended than 
the continuum emitting region. To characterise the size of the 
Bry emitting region, we fitted geometric Gaussian and ring mo¬ 
dels to the derived continuum-corrected Bry line visibilities (i.e. 
the pure Bry visibilities). We derived a HWHM Gaussian ra¬ 
dius of 0.6+0.2 mas or -0.07+0.02 AU, and ring-fit radius of 
0.7+0.2 mas or -0.08+0.02 AU for the Bry line-emitting region 
(for the adopted distance of 119 pc). 

To obtain a more physical interpretation of our AMBER ob¬ 
servations, we employed our own line radiative transfer dis c 
wind model dWeieelt et al .1120111 iGrinin & Tambovtsevall201 ll) . 
Using this approach, we computed a model that approximately 
agrees with all interferometric observables. Figures lBTl and [2] 
show that our best disc wind model MW6 (Tabled agrees with 
(1) the Bry line profile, (2) the total visibilities across the Bry 
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Garcia Lopez, R. et al.: AMBER-HR observations of HD 163296 



Fig. 3. Bry intensity distribution maps of our best disc wind model 
MW6 (see Table[3} at several radial velocities, as indicated by the white 
labels in units of km s _1 . Emission from the disc continuum and the cen¬ 
tral star is not shown. The colours represent the intensity in logarithmic 
scale in arbitrary units. 


line, (3) the continuum-corrected Bry line visibilities, and (4) 
the small observed wavelength-differential phases and closure 
phases. Therefore, our AMBER HR spectro-interferometric ob¬ 
servations of the Bry line, along with a detailed modelling of 
the line-emitting region suggest that the Bry line-emitting re¬ 
gion in HD 163296 mainly originates from a disc wind with 
an inner launching radius of 0.02 AU and an outer launching 
radius of 0.04 AU. However, Fig. [3] shows that the entire disc 
wind emitting region extends to an outer radius of ~15R, or 
~0.16 AU. The modelled disc wind emitting region is more com¬ 
pact than the inner continuum dust rim radius reported in the 
literature (R„ m ~0.19-0.45 AU; 


Tannirkulam et al 


dust rim radius reported in the 

lenistv et al.11201 ()ilEisner et ail 


stvet 

2008a b). Our modelling does 


MM 

not exclude that some of the Bry emission is produced in the 
magnetosphere. However, in addition to the emission from the 
magnetosphere, emission from a more extended region, like the 
disc wind presented here, would be required to explain the inter¬ 
ferometric observations. 


V.P.G were supported by grant of the Presidium of RAS P41. We thank the ref¬ 
eree for his/her useful comments and suggestions, which helped to improve the 
paper. 
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Fig. B.l. Comparison of the observed and modelled pure Bry line visi¬ 
bilities of our AMBER observation of HD 163296. From top to bottom: 
wavelength dependence of flux, visibilities of the first, second, and third 
baseline. In each visibility panel: (1) the observed total visibilities (red, 
green, blue, as in Fig.jT}, (2) the observed continuum-compensated pure 
Bry line visibilities (pink), and the modelled pure Bry line visibilities 
(black; model MW6, Table[3} are shown. 


Table D.l. Magnetosphere model parameters 


Parameters 

MS6a 

T e (R,)( K) 

10000 

q 

0.4 

U * (km/s) 

10 

U(r T )/UK 

1 

r c (R„) 

2.0 (0.02 AU) 

h m (R*) 

1.5 

M acc (M 0 /yr) 

lxl0~ 7 


Appendix A: Individual set of observations 

Appendix B: Pure line visibilities 

Appendix C: Examples of computed disc wind 
models 

Appendix D: Examples of computed hybrid models: 
disc wind plus magnetosphere 

Appendix D. 1: The magnestosphere 

A full description of the model employed here can be found in 
iTambovtseva et all ( 2014 1. Here, only a brief description of the 
main model parameters is presented. 

Our model considers a compact, disc-like rotating magne¬ 
tosphere of height h m through which free-falling gas reaches the 
stellar surface at some altitude near the magnetic pole. The gas 
rotational velocity component ( u ) is described by 

u(r) = Uo(r/R,) p , (D.l) 

where Uq is the rotational velocity of the gas at the magnetic 
poles, r is the distance from the star, and p is a parameter. Finally, 
we assume a dependence of the electron temperature ( T e ) of 

T e (r) — T e (/?*)exp(—rl), (D.2) 

where rl = ((r - R„)/R*) q , T e (R „) is the temperature of the gas 
near the stellar surface, and q is a parameter. 

In our hybrid model, the magnetosphere is equivalent to a 
point source at our interferometer baselines. Thus, it accounts 
for the compact and unresolved Bry emission, whereas the disc 
wind component is responsible of the resolved Bry emission. 

Because of the spread on the measured M acc values, and the 
large uncertainties in measuring this quantity (~20%), an aver¬ 
age value of 1x10 7 M 01 yr 1 was assumed in our magnetosphere 
model. 
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Fig. A.l. All three AMBER observations of HD 163296 with spectral resolution of 12 000. The three columns show the observed intereferometric 
results from 2012 May 11 (left), 2012 May 12 (middle), and 2012 June 04 (right). Shown from top to bottom are wavelength dependence of flux, 
visibilities, wavelength-differential phases (for better visibility, the differential phases of the first and last baselines are shifted by +50° and -50°, 
respectively), and closure phase observed at projected baselines as shown in the plot. The wavelength scale at the bottom is corrected to the local 
standard of rest. The typical visibilities, differential, and closure phases errors are ±5%, 5°, and 15°. 
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velocity (LSR) [km/s] velocity (LSR) [km/s] 




2.160 2.162 2.164 2.166 2.168 2.170 2.172 2.160 2.162 2.164 2.166 2.168 2.170 2.172 

X [pm] X [pm] 

Fig. C.l. Same as Fig. lB. lK left panel) and Fig.[2](right panel) but for model MW26 (see, Table[3]i. 
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velocity (LSR) [km/s] velocity (LSR) [km/s] 




2.160 2.162 2.164 2.166 2.168 2.170 2.172 2.160 2.162 2.164 2.166 2.168 2.170 2.172 

X [pm] X [pm] 

Fig. C.2. Same as Fig. lB.lI Qeft panel) and Fig. [2] (right panel) but for model MW47 (see, Table[3]i. 
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velocity (LSR) [km/s] velocity (LSR) [km/s] 




2.160 2.162 2.164 2.166 2.168 2.170 2.172 2.160 2.162 2.164 2.166 2.168 2.170 2.172 

X [pm] X [nm] 

Fig. D.l. Same as Fie. lB. ll tleft panel) and Fig.[2](right panel) but for the hybrid model MW6+MS6a (see Tablesl3land lD. It . 




































































































































































































































































